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Abstract 

We study the quantum transition at T = in the spin-^ Ising spin-glass in a 

transverse field in two dimensions. The world line path integral representation 

of this model corresponds to an effective classical system in (2+1) dimensions, 

which we study by Monte Carlo simulations. Values of the critical exponents 

are estimated by a finite-size scaling analysis. We find that the dynamical 

exponent, z, and the correlation length exponent, v, are given hy z = 1.5±0.05 

and v = 1.0 ± 0.1. Both the linear and non-linear susceptibility are found to 

diverge at the critical point. 
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Much attention has been given to the finite temperature transition in spin glass systems, 
see e.g. and reasonable agreement between theory and experiment has been obtained. 
This transition is driven by thermal fluctuations controlled by the temperature. However, 
one can also control the strength of quantum fluctuations by altering parameters in the 
system. Turning up the quantum fluctuations will decrease the transition temperature Tc, 
eventually forcing it to zero. Critical fluctuations near the transition are classical as long 
as Tc > 0, because they occur at a frequency u satisfying hu <ti ksT p. Consequently, the 
universality class is that of the classical problem except if one tunes through the transition 
at T = 0. This quantum universality class has not been much studied for the spin glass 
problem, though other quantum phase transitions, such as the metal-insulator ||^ and bose- 
glass 1^ transitions, have attracted a lot of attention. Most theoretical work on the quantum 
spin glass has been conflned to the inflnite range model, which is expected to describe 
the transition in a short range system of sufficiently high space dimension. 

Recently, however, the quantum spin glass transition was studied experimentally in an 
Ising system with dipolar couplings in which was driven to zero by applying an effective 
transverse field. Interestingly, the non-linear susceptibility, Xnu which diverges at the finite- 
T classical transition [|I|], was found not to diverge, or at least to diverge much less strongly 
than in the classical case. Furthermore, the phase transition in a quantum Ising spin system 
in (1+1) dimensions has recently been studied in detail P, see also |0. It is found that both 
the linear and non-linear susceptibility diverge not only at the critical point but also in part 
of the disordered phase. Although this model does not have frustration, and therefore might 
miss some of the spin glass physics, it is interesting to investigate whether similar behavior 
also occurs higher dimensions. It is therefore an appropriate time to study the quantum 
Ising spin glass and here we report on results of Monte Carlo simulations on a short range 
model in (2+1) dimensions. Similar calculations and analysis have also been performed in 



(3+1) dimensions pO | 



The model system studied in this paper, which is appropriate for the experimental sys- 
tem, LiHoi: Yi_^F4 , is the Ising spin glass in a transverse field with Hamiltonian 



^ = -E'^^.<^I-rE<' (1) 

where the cxj are Pauh spin matrices, F is the strength of the transverse field and the nearest 
neighbor interactions, Jij, are independent random variables with a Gaussian distribution 
of mean zero and standard deviation unity. 

If r = the Hamiltonian in is the classical two dimensional Ising spin-glass. The 
ground state is doubly degenerate (the two states being related by global spin-fiip symmetry) 
so, at T = 0, the Edwards-Anderson order parameter ^ qea = [(cf )^]av is unity. We denote 
a statistical mechanics average by angular brackets, (■■■), and an average over the quenched 
disorder by square brackets, [■ ■ -Jav Switching on the transverse field mixes the eigenstates 
of and thus diminishes the EA-order parameter, causing it to vanish at some finite value, 
Tc- This is the transition that we study here. Details of the calculations will be given 



elsewhere 12 



It is well known [T^ that the ground state energy of the (i-dimensional quantum me- 
chanical model (|1|) is equal to the free energy of a l)-dimensional classical model, where 
the extra dimension corresponds to imaginary time, i.e. 

E(T = 0) , T ^ 1 1 ^ . 
—; — - = lim — jTr e"^^ = ^Tr e""^ (2) 

where the imaginary time direction has been divided into Lr time slices of width Ar ( ArL,- = 
and the effective classical action, S, is given by 

s = -T.T.k^s,{^)s,{t) -Y.T.ks.{t)s,{t + 1) , (3) 

T {ij) T i 

where the >S'j(r) = ±1 are classical Ising spins, the indices i and j run over the sites of 
the original (i-dimensional lattice and r = 1,2,...,Lt- denotes a time slice. In Eq. (^, 
Kij = ArJij and exp(— 27^) = tanh(Arr). Note that we have the same random interactions 
in each time slice. We should take the limit Ar 0, which implies K^j and K oo. 
This extremely anisotropic limit is inconvenient for calculations but universal properties are 
expected to be independent of Ar so we take Ar = 1 and set the standard deviation of the 
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Kij to equal K. Thus K, which physically sets the relative strength of the transverse field 
and exchange terms in (|^), is like an inverse "temperature" for the effective classical model 
in®. 

We study the model (0) in c? = 2 dimensions by Monte-Carlo simulations on a simple 
cubic lattice of size L x L x L^- using periodic boundary conditions. Since various quantities 
of interest show a very strong dependence on the disorder realization we have to average 
over a large number of samples — we took 2560 samples for each temperature and size. 
The largest systems were 20 x 20 x 50, where we used up to 10^ Monte Carlo sweeps for 
equilibration plus 10^ sweeps for measurements, which were performed every 20 sweeps. 
Equilibration was checked with standard methods The simulations were performed on 
a large transputer array (GCell024 from Parsytec). 

Because the system in ® is very anisotropic, it is expected to have two different diverging 
scales: one is the correlation length in the space direction, ^ ~ 6''^, where 6 = Kc/K — 1 
is the distance from the critical point Kf.., and the other is the correlation time, in the 
(imaginary) time direction, where ~ with z the dynamical exponent. According to a 
finite size scaling hypothesis extended to anisotropic systems various thermodynamic 
quantities close to the critical point depend on two independent scaling variables, which we 
can take to be 5L^/'^ and the aspect ratio Lr/L^. The scaling analysis is straightforward 
only if it depends on a single parameter, so it is necessary to fix the aspect ratio. Since z 
is unknown, one has to scan several different sample shapes to see which choice for z scales 



best, and we follow an efficient method of doing this suggested by Huse . 

As in standard spin-glass theory we define the overlap between the configurations of 
two replicas, 1 and 2, with the same disorder as 

Q = j^T.S\'\r)S?\T) , (4) 

and for each disorder realization we calculate the dimensionless combination of moments 

<7 = 0.5[3-(Q^)/(gyl . (5) 
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The disorder averaged quantity, (^av = [s'lav [0, obeys the finite size scahng form 



g^,{K, L, L^) = ~g^,{ 5L^'\ L,/U ) , (6) 



and has the property |TT[] that it vanishes in the disordered phase for L oo, and tends to 
a finite value in the ordered phase. Consequently, g{x, y) vanishes at fixed x both for y ^ 
(where the system is a classical two-dimensional spin glass at finite "temperature", which 
is disordered) as well as for y ^ oo (where the system is effectively a long one-dimensional 
chain along the r direction, which is also disordered). Hence, g{x,y) must have a maximum 
at some value of y for fixed x. The value of this maximum decreases with increasing L in 
the disordered phase K < (where 6 = {Kc/K — 1) > 0) and increases with increasing 
L in the ordered phase. We use this criterion to estimate the critical coupling which we 
find is given by K^^ = 3.275 ± 0.025. The data are shown in Fig. 1. Furthermore, at the 
critical point, the values of L and Lr for which ^ is a maximum are related by Lr ^ L'^ . By 
this method we determine the dynamical exponent and get z = 1.50 ± 0.05. The finite-size 
scaling hypothesis can be checked a posteriori by a scaling plot for g^^ at Kc as shown 
in Fig. 2. 

Systems with fixed aspect ratio, L^/L^, can be used them to determine critical exponents 
via the usual one-parameter finite-size scaling. First of all, from Eq. (P) the derivative of 
g with respect to K at gives u and we find u = 1.0 ± 0.1, see Fig. 3. The rigorous 
inequahty u >2/d [|17| is therefore satisfied, perhaps as an equality. 

There are various susceptibilities that one can define for this problem, with different 
numbers of integrations over imaginary time. For example, the second moment of Q, Xq = 
L'^L^[((5^)]av, has a single integral over r. Defining the exponent 7q by xq ~ 5"'^'^ , then, at 
the critical point, the size dependence is given by xq ~ L?'~'^ where 7q = {2 — ifjv. On the 
other hand, the equal time spin glass correlation function, Cq = Z]i[('S'to(ro)S'j(ro))^]av, has 
no r sum and so varies as Consider next the overlap 
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l''' = JJj2i: st\n)si'\T,) . (7) 



which involves a double sum over r. The corresponding susceptibility, Xg = 
L'^L'^ [{ {q^'^y )]av; involves two time integrals so it should vary as L^^''^^ at criti- 
cality ||18[. The experimentally measured non-linear susceptibility is the fourth deriva- 
tive of the free energy with respect to a field coupling to S^, and so is related to the 
fourth order cumulant of the total magnetization by standard linear response theory, 
Xni = [{M^) — 3{M^)'^]^^ / (L'^Lr), where M = J2i,TSi{T). Since the disorder average gives 
zero unless each spin occurs an even number of times, Xni can be expressed (neglecting a 
local piece which diverges less strongly) as 



d r3 



Xni = L'^L: 



(8) 



which has three sums over r and so should diverge at criticality like "^^"^^ [|T8|. Fig. 4 



shows data and fits for Cq, Xq, Xq Xni at criticality. All the data are consistent with the 
exponent values, r] ~ 0.5, z ^ 1.5. In particular, Xni ~ L'^''^ at criticality, or equivalently 
Xni ~ -^T^ using z = 1.5. Since oc /5, Xni varies as T^^-^ for T ^ at the critical 
transverse field T^, which is quite a strong divergence. Note that, by contrast, the equal 
time correlation function does not diverge (or only does so marginally). This is because 
spatial correlations fall off quite rapidly at criticality, like r~^, as we have verified directly. 

According to scaling theory @], the (unsquared) on-site correlation function at the critical 
point, C{t) = [(5'i(0)S'j(r))]av varies as 7--{'i+2-2+»?)/{2^)^ ^j. ^-2/3 ^gj^g q^^. values for the 



exponents. Integrating this over r to get the uniform susceptibility, xf [0, Xf = Z]rC'(T), 
one finds a divergence of the form Ly^, or xf ~ T~^/^ as T ^ at F = Fc. Thus, in 
contrast to the classical spin glass [jl|, the uniform susceptibility diverges at the quantum 
spin glass transition in (2+1) dimensions. 



Similar calculations and analysis have been performed on a (3+1) dimensional model [jlO 
with results which are quite similar to ours, though the numerical values for exponents 
are somewhat different as expected. The main qualitative difference is that the uniform 
susceptibility does not diverge in (3+1) dimensions. Both our work and the results in (3+1) 
dimensions [0 show a substantial divergence of Xni, which appears to be rather different 



from experiment 0. The reason for this discrepancy is unclear at present. For future work it 
will be interesting to investigate whether the uniform and spin glass susceptibilities diverge 
in part of the disordered phase, as happens in d = (1+1) because of Griffiths singularities 
arising from rare regions which are more strongly coupled than the average ||^. 
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FIGURES 

FIG. 1. The averaged cumulant gav{K, L, Lr) for three different couphng constants (i^^^=3.20 
left, i^~^=3.30 middle and K~^=2>AQ right) and various systems sizes (L=4 (o), L=6 (+), L=8 
(□), L=12 (x) and L=16 (A)) as a function of L^- The maximum increases with L for i^~^=3.20, 
which implies K^^ > 3.20, and it decreases with increasing L for K^^ = 3.40, so K^^ < 3.40. We 
have also data for i^^^=3.25, from which we conclude that is between 3.25 and 3.30. The 
errorbars are smaller than the symbols. 

FIG. 2. A scaling plot of gi,^{K, L, L^) at K''^ = 3.30 ~ K'^ clS cl function of the scaled system 
size in the (imaginary) time direction L-t/U^^^. For each lattice size, L™^^ is chosen so that all 
the data collapses on to a single curve. The sizes are (L=4 (o), L=6(+), L=8 (□), L=12 (x) and 
L=16 (A)). The inset shows the dependence of L^^^ as a function of L. From Eq. (^) the slope is 
equal to the dynamical exponent z and a fit gives z = 1.50 it 0.05. 

FIG. 3. The derivative of Qav with respect to at = 3.30 ~ K^"^, for systems of size 
4 X 4 X 4, 6 X 6 X 8, 8 X 8 X 14, 12 x 12 x 24 and 16 x 16 x 34, which have a roughly constant 
aspect ratio, Lt/L^ , since z ~ 1.5. A least squares fit of the data to a straight line yields a slope 
of l/u = 1.0 ±0.1. 

FIG. 4. The equal time correlation function, Co, and the susceptibilites XQi Xq ^^id Xni as a 
function of L close to the critical point, = 3.30 ~ K^^, on a double logarithmic plot. The 
slopes are expected to be 2 — — z, 2 — 7], 2 — 7] + z, and 2 — r] + 2z, respectively. A least squares fit 
gives the values 0.2±0.1, 1.4ib0.1, 3.1 ±0.1 and 4.7±0.2, which are consistent with the exponents, 
r] ~ 0.5, z ~ 1.5. The system sizes are are the same as in Fig. 3. 
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Fig. 4 



